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The effects of purely elastic collisions on the dynamics of heavy inertial particles is inves- 
tigated in a three-dimensional random incompressible flow. It is shown that the statistical 
properties of inter-particle separations and relative velocities are strongly influenced by 
the occurrence of sticky elastic collisions — particle pairs undergo a large number of colli- 
sions against each other during a small time interval over which, hence, they remain close 
to each other. A theoretical framework is provided for describing and quantifying this 
phenomenon and it is substantiated by numerical simulations. Furthermore, the impact 
of hydrodynamic interactions is discussed for such a system of colliding particles. 



1. Introduction 

Dust, droplets, bubbles, and other finite-size particles suspended in turbulent flows are 
common in nature (see, e.g., Csanady 1980 Post & Abraham 2002 Shaw 2003). Their 



statistical properties are very different from that of tracers, i.e. point-like particles with 
the same mass density as the advecting fluid. Indeed, when the suspended particles 
have a finite size and a density different from that of the fluid, inertial effects become 
important. Consequently, the motion of particles starts differing from the underlying flow. 
This results in intricate correlations between the particle positions and the geometry of 
the turbulent flow. Heavier particles are expelled from vortical structures while lighter 
particles concentrate in their cores. A consequence of these mechanisms is the presence 
of strong fluctuations in the spatial distribution of particles. This phenomenon, known as 
preferential concentration, has been the subject of extensive research in fluid dynamics 
for the last decades (see Douady et aL|1991 Squires fc Eaton|[l991 ). 

Another consequence of inertia is that particles are likely to be very close to each 
other with large velocity differences. The process leading to such events, known either as 



2005 


I, has been 


Bee et al. 


2010 



e-g-, 



is known to strongly alter possible interactions between particles. Since the pioneering 



work of Saffman & Turner ( 1956 ) motivated by coalescences of cloud droplets, much work 
has been devoted to understanding the rate at which heavy inertial particles collide. A 
commonly adopted approach consists in assuming that the inter-collision time is much 
longer than the convergence timescale of particle dynamics to a statistically stationary 
regime. This premise, which is asymptotically true in the limit of very dilute suspensions, 
entitles counting collisions without having to effectively perform them. The frequency at 
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which such ghost particles collide is then a time-independent statistical observables that 
can be quantified as a function of particle sizes and response times. Several studies have 
assumed this ghost-particle hypothesis in order to estimate collision rates of heavy inertial 



particles in turbulent flows (see, e.g., Sundaram & Collins 1997 Falkovich et al. 2002) 



However, little is known about the limits of such an approach. For instance, an important 
statistical weight is given to events when two concentrated clouds of particles cross each 
other with a large velocity difference. It is clear that the very-dilute approximation should 
then fail and that multiple-collisions are likely to occur. We report in this paper results 
on the effects of actual collisions on the dynamics and statistics of inertial particles 
transported by a non-stationary fluid flow. 

The simplest framework for treating short-range collective effects in the particle dy- 
namics is to consider an ensemble of hard spheres that are suspended in a prescribed flow 
and which undergo purely elastic (momentum and energy preserving) collisions with each 
other. Such a system, in the absence of any underlying fluid transport, has long been 
a paradigmatic model in statistical mechanics. It is the basis of the kinetic theory of 
gasses, in which it is assumed that the kinetic energy is conserved both by the collisions 
and by the dynamics of each particle. If we allow for a certain amount of inelasticity in 
the collisions a suitable model for granular gasses is obtained (see Goldhirsch 2003, for 
a review). A peculiar behaviour of such granular systems is the spontaneous aggregation 
of particles. Since inelasticity implies the loss of a finite percentage of kinetic energy at 
each collision, the particles can eventually stick together leading to the formation of large 
clusters. However, this effect occurs also in settings that are not common in studies of 
granular media. Indeed, as we will see in this work, it suffices that kinetic energy dis- 
sipation occurs not via collisions but rather through individual particle dynamics. This 
is the case for heavy inertial particles whose motion is dominated by viscous damping 
and which undergo purely elastic collisions. We show that the clustering phenomenon 
emerging in such systems originates from what we call sticky elastic collisions. During 
these events, the particles bounce many times against each other and energy is dissipated 
during their motion between successive collisions. This mechanism has strong influences 
on the statistics of inter-particle distances and relative velocities. Furthermore, to vali- 
date the presence of this effect in real settings, we study the influence of hydrodynamic 
interactions and show that they cannot prevent sticky elastic collisions from occurring. 

This paper is organised as follows. In we discuss the equations of motion of the iner- 
tial particles and recall some key results, in the absence of any collisions, on preferential 
concentration. In fj3| we present numerical results and provide an asymptotic analysis for 
the phenomenon of sticky elastic collisions in the absence of hydrodynamic interactions. 
We then discuss the effect of hydrodynamic interactions in fj2] by considering far-field 
interactions. We make some concluding remarks and summarise our results, as well as 
provide a perspective for future work, in S|5| 



2. The Model 

To set the stage, we begin by recalling the basic physics of a system of N small hard 
spheres which are in a random, time-dependent, incompressible fluid field u(x, t) and are 
subject to viscous dissipation. When such particles have a small Reynolds number and 
are much heavier than the fluid, they interact with the flow by Stokes viscous drag and 
their trajectories Xi(t) are determined by Newton's law: 



Xi = Vi, Vi = 



— [ Vi -u(xi,t)] ie[l,N] 



(2.1) 
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where r is the viscous-drag relaxation (Stokes) time defined via r = 2p p a? / (Qpfv), where 
p p is the particle density, p f the fluid density, v its kinematic viscosity, and a is the particle 
radius. The Stokes number is a measure of the inertia of the particle and is defined as 
St = t/t{, where Tf is a characteristic timescale of the fluid flow. In this work, we 
additionally introduce interactions between particles in the following sense. The particles 
interact through elastic collisions which, for the case of spherical particles of equal size 
and mass, correspond to an exchange of the radial component of the velocities of the two 
colliding particles upon impact, namely when | 

The clustering phenomenon occurs naturally, independently of the carrier flow com- 



pressibility, for an ensemble of particles which evolve according to Eq. (2.1 1, even in the 
absence of interactions or collisions between particles. We note that the case of no colli- 
sions is equivalent to the ghost particle approach. The physical mechanisms which lead to 
strong inhomogeneities in the spatial distribution of particles arise from the underlying 
dissipative chaotic dynamics: The system is characterised by a constant contraction rate 
d/r in the position- velocity phase space which drives the particles towards a dynamically 
evolving fractal set The clusters of particles are the projection of such fractals on the 
position space. The fractal dimensions of particles distribution provide a convenient tool 
to quantify their clustering (see Bee et al. 2007 Calzavarini et al. 2008[ ). In particular 



the correlation dimension Di is defined via the power-law behaviour of the probability 
distribution function (PDF) of inter-particle distances Vnif) ~ r D2_1 (see the inset of 
Fig.[l^i for particles with different Stokes numbers in a three-dimensional random flow). 
Note that the density P2 (r) is such that the probability that two particles are at a dis- 
tance between r and r + dr is P2(r) dr and relates to the radial distribution function g(r): 
in three dimensions, one has g(r) = 4irr 2 P2(r). 

The intensity of clustering is influenced by the properties of the velocity field u(x, t) 
and in particular by its spatial and temporal correlations. Here we will assume that the 
velocity is differentiable in space and time and characterised by unique time and length 
scales. Let us denote by Lf the fluid flow correlation length, which is assumed much larger 
than the particle radius a, and by Tf its correlation time, which is of the same order as 
the turnover time L[/U where U is the typical amplitude of u. In the limit of vanishing 



dissipation St — t/ti — > 0, Eq. (2.1) becomes that of tracers, namely Xi — u(xi,t), 
and the incompressibility condition V ■ u — ensures a uniform distribution of particles. 
Particles distribute uniformly also in the opposite limit St — > oo, in which the force 
acting on particles become very small so that they follow almost a ballistic motion and 
fill the whole position-velocity phase-space. The maximum of clustering is achieved for 
intermediate (order-unity) values of St where the minimum of the fractal dimension is 
reached. 

Clearly, the presence of collisions will affect the particle distribution on scales of the 
order of their size. Heuristically, one expects the two-particle distribution P2(r) to be 
unchanged at separations r much larger than the interaction distance 2 a. At separations 
of the order of 2 a, the elastic collisions decorrelate the particle dynamics from the carrier 
flow. Naively, one would then expect that the particles distribute like an ideal gas at such 
scales, and consequently have a uniform distribution. In dilute systems, the crossover 
between these two regimes would occur at a scale given by the distance travelled by the 
particles before it relaxes to its attractor. This distance can be written as r* = Tv c (St), 
where v c is the typical velocity difference at collisions. However, as seen from the inset of 
Fig.[l^, this naive picture seems wrong. As we will see in the next section, we indeed find 
that the two-particle density p2 diverges when r — > 2 a. This effect is due to the presence 
of sticky elastic collisions. 
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1. (a) Log-log plot of the cumulative probability distribution P^ir) as a function of 
for some representative values of St as indicated in the legend. Inset: Log- log plot of 
the two-particle density P2(r) = dP^(r)/dr as a function of r for the same values of St. (b) 
Exponent a versus St. The gray dashed line is an empirical fit of the form — exp( — (St/0.09) ). 
Inset: local slopes qls = d\ogP^(r — 2a)/dlog(r — 2a) — 1 as a function of (r — 2a) for various 
St. The two dashed vertical lines indicate the region over which we calculate the mean and the 
standard deviation of q:ls- 



3. Particle adhesion through recurrent collisions 

To investigate the effect of collisions on clustering we resort to numerical simulations 
of Eq. (2.1 ). For simplicity and without any loss of generality, we consider two (N = 2) 
particles in a three-dimensional cubic domain of size L with periodic boundary conditions. 
The velocity field is obtained from a superposition of Fourier modes whose amplitudes 
are stochastic Ornstein-Uhlenbeck processes with Gaussian statistics and a correlation 
time Tf. The amplitude of each mode, which has a standard deviation of the order of 
L/t{, is chosen to ensure statistical isotropy at small scales. In our simulations we use 
several values of r with Stokes numbers St = t/t{ lying between 0.01 and 1.1. We have 
done simulations with various values of the particle radius a, but report here results 
obtained for a 0.02. Our time marching is an implicit Euler scheme with a fixed time 
step St = 10 -4 when the particles are far away from each other. However, when the 
particle are close to each other, the time step is adapted in order to resolve collisions 
with a high accuracy. 

Our numerical simulations reveal a very interesting, hitherto unknown, phenomenon 
which is clearly absent in the collisionless or ghost-collision case. As expected the effects 
of collisions on the spatial distribution of particles, as characterised by the PDF of inter- 
particle distance, are negligible for particle separations much larger than 2 a. However 
collisions dramatically affect the statistics of pair-separations at small scales. In partic- 
ular, the probability distribution function of the inter-particle distance displays a 
power law behaviour P2 (r) ~ (r — 2a) a for distances close to the cutoff 2 a. The exponent 
a is a monotonically increasing function in the Stokes number St : It begins with the 
value — 1 for St — > and approaches as St — > 00 . Figure [T^i shows the cumulative 
probability distribution P^ir), which is the probability of finding two particles at a dis- 
tance less than r, as a function of (r — 2a), on a log-log scale, for some representative 
values of St that we use in our simulations. We clearly find that P 2 < behaves as a power 
law oc (r — 2a) a+1 at small values of r — 2 a. From a local slope analysis we extract 
the local scaling exponent ols = d log P 2 < ( r ~ 2a)/dlog(r — 2a) — 1 (see the inset of 
Fig. [l]}); the mean of this gives us a measure of the scaling exponent a and the standard 
deviation an estimate of the error. In Fig. [lb we show the behaviour of a as a function 
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of the Stokes number. The gray dashed line is a fit to our numerical data and is given by 
— exp(— (5t/0.09) 2 ). In the inset the dashed vertical lines denote the region over which 
we calculate the mean and the standard deviation to obtain a. Two asymptotic values 
are clearly visible on Fig. [TJd: one observes that a — > when St — > oo and a — > — 1 
when St —> 0. Note that the stiffness of the system in the limit St — > prevents us 
from obtaining accurate numerical results in this limit. Nevertheless, our data show a 
monotonic convergence towards a = — I. These two asymptotic values of the exponent 
are signatures of two different collision mechanisms. 

For large Stokes numbers the motion of particles is weakly correlated with the local 
values of the fluid velocity field, because it is determined by the cumulative contributions 
of the flow integrated over the particle trajectories with a long memory kernel. Close 
to the collision distance individual particle velocities arc almost uncorrelatcd and vary 
on very large timescales. The relative motion of the two particles is therefore almost 
ballistic, with a random relative velocity. The time that the two particles spend at a 
distance between r and r + dr is given by dt = V dr, where V is a their typical velocity 
difference. This leads to p%{r) ~ (r — 2 a) thus yielding a — > when St — > oo. 

The limit St — > is more complicated. In this case the motion of the particles is 
strongly correlated with the fluid velocity field u, and therefore particles typically arrive 
at colliding distances with a very small relative velocity, of the order of the fluid velocity 
difference, namely 2 a a, where a is the local gradient of u. After the collision they tend 
to separate but the fluid velocity field quickly brings them back together, because of 
the short relaxation time r < Tf. Particles then collide again and this mechanism will, 
therefore, lead to a long series of high-frequency collisions during which the particle 
remain within a distance of the order of 2 a. Next, to quantify the effects of these events 
on the statistics of inter-particle distance, let us consider a simple, one dimensional model 
for the separation r = \x% — and the radial relative velocity v — (vi — V2) ■ f of two 
particles, between two collisions occurring at time t n and t n+ \. The equations of motion 
can be written as 

r = v , tv = —v — 2ao. (3-1) 

Here, we have assumed that the fluid velocity gradient —a < remains constant. This is 
justified because particle dynamics is much faster than the correlation time of the fluid 
velocity (r <C Tf). Also, we have neglected higher-order terms in the Taylor expansion 
of u. The above equation gives r(t) — 2a — t(v u + 2aa)[exp(— t/r) — I] — 2aat and v(t) = 
v n exp(— t/r) + 2acr[exp(— i/-r) — 1] where v n is the relative radial velocity immediately 
after the n-th impact at time t n . Introducing the small parameter e„ = v n /(2ao~) and 
Taylor expanding the above expressions for f/r < 1 one obtains a recursive relation 
for the relative velocity at collision, namely e n +i = e n (l ~ 2e n /3), which leads to v n ~ 
2aa/n. The inter-collision time 6 n = t n+ i —t n decreases as 8 n ~ t/ti, and the maximum 
distance r* reached by the two particles in the excursion between two collisions scales 
as 5 = r*/(2a) — I ~ ar/n 2 . The number of collisions increases exponentially in time 
n c ~ exp(Ci/r), with C > 0, and the inter-particles distance goes to 0. We call this 
phenomenon sticky elastic collisions. Note that the series defined by the sum of inter- 
collision times J2 n ®n is not converging. This ensures that the number of collisions does 
not become infinite in a finite time, at variance with the case for inelastic collisions 
where the particles eventually collapse and aggregate. In the case of the one-dimensional 
sticky elastic collisions, the recurrent process stops only when the fluid velocity gradient 
becomes positive (for t ~ Tf) and takes the two particles far away. During any one of 
these events, only the first m collisions will contribute to the probability of having the 
two particles at a distance larger than r ~ 2a(l + aT/m 2 ). The fraction of time spent at a 
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Figure 2. (a) Log-log plot of the PDF of inter-collision radial velocity v c for St = 0.01; the 
dashed line shows a scaling of —1.85. (b) Log-log plot of the PDF of inter-collision time 6 for 
St — 0.01; the dashed line shows a scaling of —1.95. 



distance larger that r is thus ~ J2 n <m 9 n ^ \nm ^ — ln(r — 2a). This gives for the PDF 
of inter-particle distances P2(r) ~ (r — 2a) _1 , yielding the asymptotic value a — > —1 for 
St -> 0. 

The heuristic arguments developed here to quantify the statistical signature of sticky 
events are purely one-dimensional. However, they extend to higher dimensions by geomet- 
rical considerations. The one-dimensional case would hold true if the velocity difference 
v between the particles were exactly aligned with their separation r. A misalignment 
leads to rebounds of the two particles at different locations on their surfaces. For spher- 
ical particles, this implies that better is the alignment between v and r, higher is the 
number of successive secondary collisions. Statistically, this implies that the distribu- 
tion of distances is dominated by almost head-on collisions, which can essentially be 
treated as a one-dimensional problem. The space dimensionality should just appear as a 
multiplicative factor in the power-law behaviour at r — > 2 a. 

The effects of these events are detectable also in the statistics of inter-collision times 9 
and relative radial velocity at collision v c . Within the model derived above, for St — > 0, 
we obtain p(v c ) ~ v~ 2 and p{9) ~ 9~ 2 . In Fig. and Fig.fJ^) we show a log- log plots of 
p(v c ) versus v c and p(0) versus 9, respectively, obtained from simulations for St = 0.01; 
the dashed lines shows a scaling exponent of ss —1.9 in each case. We have checked 
the scaling becomes significantly shallower and the extent of scaling gets progressively 
reduced with St increases. Note that a similarly steep increase of the inter-collision time 
distribution at small values has been observed by Ten Cate et al. (2004) in full direct 
numerical simulations of finite-size particles suspended in a viscous flow. This effect has 
been interpreted as a possible consequence of lubrication forces between the particles that 
makes the particles remain close to each other for long times. In the next section, we 
will comment on the influence of hydrodynamical interactions on sticky elastic collisions. 
Nevertheless it is worth stressing here that lubrication is not necessary to obtain multiple 
collisions between particles. 

The tails of p(v c ) and p(9) at small values seem to indicate that they cannot be 
normalised. However, we observe for small but finite values of St, a cutoff at the smallest 
values of v c and 9, which prevents this divergence. Indeed, the power-law behaviours are 
due to typical sticky events. As we have seen, the number of collisions is of the order 
of n c ~ exp(C/St). Therefore, the minimal collisional velocity and inter-collision time 
are both oc l/n c ~ exp(—C/St). The power-law is thus just appearing in intermediate 
ranges, namely (2 a/n) exp(-C/St) < v c < (2o/r f ) and rexp(-C/5i) < 9 < r. The 
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events leading to values smaller than the lower bounds are related to situations where 
the fluid velocity gradient is maintained negative for an exceptionally long time. This 
happens with an inverse-Gaussian probability in a finite-correlation-time flow, whence 
the cutoff. Also, the extension of this argument to the singular limit St = is far from 
obvious. Note finally that in dimensions higher than one, the geometrical considerations 
explained above should also provide a cut-off for very small collision times and velocities. 



4. Effect of hydrodynamic interactions 

In this section, we address the question whether or not sticky elastic collisions can be 
observed in realistic flows. For that, it is important to know how far the observations and 
conclusions drawn above are valid when we, in addition to collisions, introduce hydro- 
dynamic interactions between particles. To answer this question, we take into account 
the far-field or long-range interactions and assume that it is valid all the way to the 
smallest separations. This approach yields several interesting results which we discuss 
below. Of course this assumption, in reality, breaks down when particles approach very 
close to one other. However in that case, either lubrication leads to an effective increase 
in the particle physical radius or, in the case when the velocity difference is too large, 
hydrodynamics might not be a valid description of the interactions between the particles. 

In its simplest formulation, the long-range hydrodynamic interactions between parti- 
cles in a flow is taken into account by considering the perturbation in the ambient fluid 
velocity field, as experienced by an individual particle, because of the motion of all the 
other particles. Thus the effective velocity field acting on any particle is a superposition 
of the unperturbed (turbulent) advecting flow u(x, t) and of the perturbation to this flow 
due to the other particles. Formally, this perturbation u^' on the isolated ith particle 
due to another particle j of radius a at a distance \r^ \ and moving with a velocity 
can be written as a combination of a Stokeslet and a potential dipole flow as follows (see 
Wang et aL|[2009| : 



3 a 3 a 3 
4^0 ~ 4 |r(j)| 3 



| r (j)| 2 



3 a 1 



4|r0')| 4|r(j)| 3 



v&>. (4.1) 



In a system of N particles, the net perturbation on the flow field experienced by any 
particle i is obtained by summing over the contributions made by each of the other 
(N — 1) particles. Given the structure of the equations, for a system of N particles, it 
is impossible to solve exactly the perturbation field. Thus various approximations and 
iterative schemes become essential. However, in the present problem being studied in this 
paper, which involves two particles only, it is possible to solve exactly the hydrodynamic 
interaction term as it involves merely an inversion of a 6 x 6 matrix. 

Let us now try to understand the effect of hydrodynamic interactions on the particle 
dynamics from a theoretical point of view. Without any loss of generality, and considering 
only two-particle interactions, let us consider a model where the two particles approach 
each other with the same velocity v. A direct consequence of this model is that the 
perturbation on particle 1 due to particle 2 is equal and opposite to the perturbation 
it' 2 ) on particle 2 due to particle 1, i.e. = —u^ 2 \ The equation of motion for particle 1 
can be written as 



~d7 ' r 



- (v - m - tt {2) ) ; (4.2) 



1 m\ /3 a la 3 



„(») = ( w + i„-«W)(~-~J l (4.8) 
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Figure 3. (a) Exponents a of the two-particle distribution pz(r) oc (r-2o)", along with their 
error bars. The case with hydrodynamic interactions (x) is represented versus St e ff — (5/16) St 
and that without interactions (•), as a function of St; the gray dashed line is the empirical fit 
-exp(-(5"i/0.09) 2 ) discussed in previous section, (b) PDF of inter-collision radial velocity v c ob- 
tained from simulations with hydrodynamic interactions (x) for St = 0.15, that is St e a = 0.047 
and without hydrodynamic interactions (•) for St = 0.05. 



where er denotes the unperturbed fluid velocity gradient. Since, at the point of collision, 
\r\ = 2a, one has 

=-— (v- 2cr fa). (4.4) 
5 



By using Eq. (4.4) in Eq. (4.3), we eventually obtain 

dv 16 . , , 

_ = -_(„-2<r*«). (4.5) 

The above analysis shows that the effect of long-range hydrodynamic interactions re- 
duces in the vicinity of collisions to the dynamics in absence of interactions but with 
an effective Stokes number, which is equal to the actual Stokes number reduced by a 
factor of 16/5 = 3.2. Thus a system of particles with a Stokes number St and subject 
to hydrodynamic interactions can be replaced by a system of particles, without any hy- 
drodynamic interactions but with an effective Stokes number St c s = (5/16) St when we 
consider their statistical properties for very small inter-particle separations. 

To confirm the arguments presented above we resort once more to numerical simu- 



lations by implementing the long-range hydrodynamic interactions (4.1). We use values 
of St between 0.05 and 1.1 as we had used for the case without any interaction terms. 
We begin by measuring the values of the exponent a introduced in the previous Section 
that describes the behaviour of the inter-particle distance distribution. Figure (|3^) shows 
the values of a obtained as a function of the effective Stokes number St e g = (5/ 16) St 
(crosses); our data seems to fall, within error bars, on the empirical fit shown by a 
dashed gray line. To make the comparison more illuminating we plot on the same graph 
and as a function of the actual Stokes number St, the values of a (black dots) obtained 
from numerical simulations without hydrodynamic interactions, and already shown in 
Fig. ^p). Furthermore, in Fig. ^jp) we show the PDF of the inter-collision velocity p(v c ) 
as a function of v c , on a log-log scale, for St — 0.05 obtained from a simulation without 
hydrodynamic interactions (black dots) and for St — 0.15 obtained from a simulation 
with hydrodynamic interactions (crosses). The two PDFs are nearly overlapping as our 
arguments before would suggest that the effective Stokes is St e g = 0.047 ~ 0.05 for the 
case with the hydrodynamic interactions. 
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Our results suggest that hydrodynamic interactions increase the efficiency of dissipa- 
tive mechanisms in terms of a reduction of the effective Stokes number. However, such 
considerations can only lead to qualitative deductions as our study account for long-range 
interactions only. The effect of lubrication forces will become dominant for particles at 
very small separations. On the one hand, this type of interaction is expected to decrease 
the collision efficiency between particles (see, e.g., Wang et aZ.|2005 1 . On the other hand, 
lubrication is expected to increase damping when particles get close to each other; this 
effect is usually modelled by a restitution coefficient less than unity. Because of these 
two competing mechanisms, it is difficult predict whether short-range hydrodynamic in- 
teractions will enhance or diminish the sticky elastic collision phenomenon. 



5. Conclusions 

In our work we have considered the effect of elastic collisions on the clustering of 
inertial particles. In particular we have investigated their influence on the probability 
distribution of inter-particle distance. Surprisingly, our findings differ markedly from 
the naive picture that collisions might only introduce a small-scale molecular chaos. 
We observe that the small-distance statistics is dominated by a phenomenon, which 
we call sticky elastic collisions, during which particles undergo a very large number of 
collisions during a time of the order of the fluid correlation time. It is interesting to note 
that these sticky elastic collisions remarkably resembles inelastic collapses observed in 
granular media, even though the underlying assumption in granular media (conservative 
inter-collision dynamics and dissipativc collisions) is exactly the opposite of what we have 
considered here. In addition we have investigated the effect on this phenomenon of long- 
range hydrodynamic interactions between particles. Our results seem to indicate that 
the most significant effect at small scales of such interactions is to introduce an effective 
Stokes number. The problem of investigating the effect on sticky elastic collisions of short- 
range hydrodynamical interactions requires more rigorous theoretical understanding and 
more elaborate numerical simulations. 

In this study we have focused on two-particle interactions. It is clear that the phe- 
nomena of sticky elastic collisions will be present even for large numbers of interacting 
particles. In this light, collective phenomena may emerge and provide a new mechanism 
to dissipate kinetic energy in violent collisions between particles. There arc still many 
open questions concerning the stability of coalescence processes for the high impact veloc- 
ities that are observed in turbulent settings. In particular, estimates on relative velocities 
between meter-sized objects in circum-stellar disks are by far too large to allow for their 

). The dissipative mech- 
anisms relating to sticky elastic collisions might play a role there. 
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accretion and growth to form planet embryos (Wurm et a?.|2001 



REFERENCES 

BEC, J., BlFERALE, L., BOFFETTA, G., CENCINI, M., LANOTTE, A.S., MUSACCHIO, S. & 

TOSCHI, F. 2007 Heavy particle concentration in turbulence at dissipative and inertial 
scales. Phys. Rev. Lett. 98, 084502. 



10 



J. Bee, S. Musacchio, and S. S. Ray 



Bec, J., Biferale, L., Cencini, M., Lanotte, A.S. & Toschi, F. 2010 Intermittency in the 
velocity distribution of heavy particles in turbulence. J. Fluid Mech. 646, 527-536. 

Calzavarini, E., Kerscher, M., Lohse, D. & Toschi, F. 2008 Dimensionality and morphol- 
ogy of particle and bubble clusters in turbulent flow. J. Fluid Mech. 607, 13-24. 

Csanady, G. 1980 Turbulent diffusion in the environment. Geophysics and Astrophysics Mono- 
graphs Vol. 3 D. Reidel Publishing Company. 

Douady, S., Couder, Y. & Brachet, M.-E. 1991 Direct observation of the intermittency of 
intense vorticity filaments in turbulence. Phys. Rev. Lett. 67, 983-986. 

Falkovich, G., Fouxon, A. & Stepanov, M. 2002 Acceleration of rain initiation by cloud 
turbulence. Nature 419, 151-154. 

Falkovich, G. & Pumir, A. 2007 Sling effect in collision of water droplet in turbulent clouds. 
J. Atm. Sci. 64, 4497-4505. 

Goldhirsch, I. 2003 Rapid granular flows. Ann. Rev. Fluid Mech. 35, 267-293. 

Jeffrey, D. J. & Onishi, Y. 1984 Calculation of the resistance and mobility functions for two 
unequal rigid spheres in low-Reynolds-number flow. J. Fluid Mech. 139, 261-290. 

POST, S. & Abraham, J. 2002 Modeling the outcome of drop-drop collisions in Diesel sprays. 
Intl. J. Multiphase Flow 28, 997-1019. 

Ten Cate, A., Derksen, J.J., Portela, L.M. & Van Den Akker, H.E.A. 2004 Fully 
resolved simulations of colliding monodisperse spheres in forced isotropic turbulence, J. 
Fluid Mech.—/ 519, 233-271. 

Saffman, P.G. & Turner, J.S. 1956 On the collision of drops in turbulent clouds. J. Fluid 
Mech. 1, 16-30. 

Salazar, J.P.L.C. and Collins, L.R. 2012 Inertial particle relative velocity statistics in 

homogeneous isotropic turbulence. J. Fluid Mech. 696, 45-??. 
Shaw, R.A. 2003 Particle-turbulence interactions in atmospheric clouds. Ann. Rev. Fluid 

Mech. 35, 183-227. 

Squires, K.D. & Eaton, J.K. 1991 Preferential concentration of particles by turbulence. Phys. 
Fluids A 3, 1169-1178. 

Sundaram, S. & Collins, L.R. 1997 Collision statistics in an isotropic particle-laden turbulent 

suspension. Part 1. Direct numerical simulations. J. Fluid Mech. 335, 75-109. 
Wang, L.P., Ayala, O., Kasprzak, S.E. & Grabowski, W.W. 2005 Theoretical formulation 

of collision rate and collision efficiency of hydrodynamically interacting cloud droplets in 

turbulent atmosphere, J. Atmos. Sci. 62, 2433-2450. 
Wang, L.-P., Rosa, B., Gao, H., He, G & Jin, G 2009 Turbulent collision of inertial particles: 

Point-particle based, hybrid simulations and beyond. Intl. J. of Multiphase Flow 35, 854- 

867. 

Wilkinson, M. & Mehlig, B. 2005 Caustics in turbulent aerosols, Europhys. Lett. 71, 186-192. 
Wurm, G., Blum, J. & Colwell, J.E., A new mechanism relevant to the formation of plan- 
etesimals in the solar nebula, Icarus 151, 318-321. 



